Magnetic field induced localization 
in a two-dimensional superconducting wire network 
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We report transport measurements on superconducting wire networks which provide the first ex- 
perimental evidence of a new localization phenomenon induced by magnetic field on a 2D periodic 
structure. In the case of a superconducting wave function this phenomenon manifests itself as a 
depression of the network critical current and of the superconducting transition temperature at a 
half magnetic flux quantum per tile. In addition, the strong broadening of the resistive transi- 
tion observed at this field is consistent with enhanced phase fluctuations due to this localization 
mechanism. 



PACS numbers: 72.15, 73.23, 74.25 

In a recent paper Q| a novel case of extreme localiza- 
tion induced by a transverse magnetic field was predicted 
for non interacting electrons in a two-dimensional (2D) 
periodic structure. This new phenomenon, due to a sub- 
tle interplay between lattice geometry and the magnetic 
field, differs from Anderson localization on two essen- 
tial points: it occurs in a pure system, without disorder, 
and the system eigenstates are not localized but non- 
dispersive states. In a tight-binding (TB) approach, it 
can be simply understood in terms of Aharonov-Bohm 
effect which, at half a flux quantum per unit tile (half- 
flux), leads to fully destructive quantum interferences. 
For this flux, the set of sites visited by an initially lo- 
calized wave-packet will be bounded in Aharonov-Bohm 
cages Q. This effect is absent on other regular periodic 
lattices at half-flux, such as the square and the triangular 
lattices. 

Superconducting wire networks are suitable to address 
phase interference phenomena driven by a magnetic field 
||. These systems are extremely sensitive to phase co- 
herence of the superconducting order parameter over the 
network sites which is exclusively determined by the com- 
petition between the external field and the network geom- 
etry Hi ■ Besides, the quantum regime is accessible even in 
low T c diffusive superconductors: since all Cooper pairs 
condense in a quantum state, the relevant wavelength is 
associated with the macroscopic superfluid velocity and 
can be much larger than the lattice elementary cell jij. 
Also, the magnetic field corresponding to one supercon- 
ducting flux quantum <I> = hc/2e, is easily accessible: 
it is about 1 mT for a network cell of 1 /im 2 , in con- 
trast to the unattainable 10 3 T for an atomic lattice. 
In addition, some features of the TB spectrum, namely 
the Hofstadter butterfly || , are experimentally accessible 



in the model system of a superconducting wire network 
H||. As shown by de Gennes and Alexander [@,||, the 
linearized Ginzburg-Landau (GL) equations for a super- 
conducting wire network can be mapped onto the eigen- 
values equation of a TB hamiltonian for the same geom- 
etry. This mapping is of particular relevance since one of 
the remarkable findings of Ref. JlJ is the total absence of 
dispersion in the TB spectrum at half-flux. In the con- 
text of a superconducting network, the localization effect 
is expressed by the inability of the superconducting wave 
function to carry phase information throughout the net- 
work and therefore, transport anomalies are expected. 

In this Letter we present transport measurements on 
2D superconducting networks with the so-called T 3 ge- 
ometry (see inset of Fig. Q). Our results allow us to con- 
firm some of the exotic features of the T3 energy spec- 
trum related to the localization mechanism. The field- 
temperature (H,T) superconducting transition line is de- 
termined and related to the ground state of the T3 spec- 
trum. We also compare the critical current as a function 
of the magnetic field with calculations of the group veloc- 
ity. The striking behavior found at half-flux is discussed 
as a possible signature of localization effects. The strong 
broadening of the normal to superconductor transition 
supports this interpretation. Very few experiments were 
reported so far on localization phenomena in supercon- 
ducting networks and only the issues of irrational mag- 
netic flux H or disorder have been addressed. 

The networks pattern was defined on a 600 nm thick 
layer of positive UV3 resist using an e-beam writer Le- 
ica VB6-HR. A 100 nm thick layer of pure aluminum 
was e-beam evaporated in a ultra-high vacuum cham- 
ber, followed by the resist lift-off We designed 
two series with a large patterned area: Star 600 de- 



1 



fined on a 0.6 x 1 mm 2 surface and Star 20 denned on 
a surface of 0.02 x 1 mm 2 , which required stitching of 
200 x 200 /im 2 writing fields. The elementary tile side 
length is a = 1 /jm, the wires having 100 nm width and 
100 nm thickness. 
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FIG. 1. Reduced critical temperature vs frustration for 
Star 600 (heavy line, left axis). The theoretical curve has 
been offset by -0.001 for clarity (small dots, right axis). In- 
set: Electron micrography of the T3 network (a = 1.0 /j,m). 
Measuring current is applied along the horizontal direction. 

Dynamic resistance measurements were performed us- 
ing a 33 Hz ac four terminal resistance bridge with an 
ac measuring current of 20 nA. Sample probes were con- 
nected to the cryostat terminals by ultrasonic bonding 
of 25 /im gold wires. Non invasive voltage probes were 
placed at 0.2 mm from the current pads. The zero field 
transition temperatures 7^(0) were 1.234 K for Star 600 
and 1.240 K for Star 20, using a resistance criteria of 
half the normal state resistance R n at 1.25 K, which are 
4.20 £1 and 63.56 fl, respectively. The resistive transition 
width in zero field is 3 mK (10% — 90%) for both samples 
indicating a good homogeneity of the networks. 

The field dependent transition temperature T C (H) was 
monitored by locking the temperature controller to keep 
the sample resistance at 0.5i? n as the magnetic field is 
varied. The experimental data is to be compared with 
the lowest energy solution of the network linear GL equa- 
tions, that is given in terms of the ground state eigenvalue 
e g (f) of the TB spectrum § by, 



1 



TJJ) = g(0) 2 
T c (0) a 2 



\V18J 



(1) 



where £(0) is the superconducting coherence length at 
zero temperature and / the frustration. Neglecting field 
screening effects, / — &/§ where $ = Hc?\fZj2 is the 
magnetic flux through a rhombus tile. 

The transition line of Star 600 is plotted in Fig. |l| in 
reduced units 1 — T c (/)/T c (0) as a function of frustra- 
tion. Since the transition line is periodic on <£> we only 
displayed it in the field range < / < 1 |. A small 



parabolic background due to field penetration in the 
wires was subtracted from the experimental T c (f). We 
also display the theoretical T c (/) obtained using Eq. (|l|) 
and £(0) = 157 nm, the only adjustable parameter. The 
fine field structure of the experimental data is very well 
described by the theoretical curve. Distinct downward 
cusps are visible at low order rationals / = 1/q, for 
q = 3,4,6, and 2/9. They reflect the long range phase 
ordering of the order parameter among network sites, es- 
tablished at fields commensurate to the underlying lat- 
tice. These features were discussed previously The 
novel feature of the transition line occurs at / = 1/2, 
where the maximum of T c (f) depression (30 mK) is 
achieved, associated with an inversion of the field modu- 
lation concavity. This anomalous cusp persists distinctly 
at all criteria used on T c (/) determination, from 0.06i?„ 
to 0.87i?„, though the downward cusps at other rationals 
fade out with increasing temperature. This cusp is sim- 
ilar to the T c variation in a single loop geometry close 
to / = 1/2 and is characteristic of quantum effects 
determined on a finite length scale. It indicates that at 
half-flux the network transition is determined by fluxoid 
quantization at independent tiles. 

The maximum depression of T c (f) at half-flux shows 
the strong incommensurability at this field. To our 
knowledge, these results are the first experimental ob- 
servation of such an effect on an extended periodic net- 
work. Besides, they indicate that 2D periodicity is not a 
sufficient condition for a commensurate state to exist at 
rational /. 




FIG. 2. Transition width vs frustration for Star 600. The 
large broadening close to / = 1/2 indicates the presence of 
strong phase fluctuations. 

We also observed a strong broadening of the resistive 
transition AT W idth, at half-flux, as displayed in Fig. for 



Star 600. AT W 
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is obtained as the difference between 



the T c (/) curves taken for criteria 0.6R n and 0.1i? n , re- 
spectively. The anomalous enhancement (up to 12 mK) 
at half-flux, twice the average width over most of the field 
range, confirms the singular behavior found at this field. 
At the strong commensurate fields / = 0, 1 
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is sharply reduced to a few rnK as expected for a phase 
ordered system. Close to these fields, the phase of the or- 
der parameter at the network sites is able to " lock" in the 
nearest commensurate state with the creation of few mo- 
bile defects, broadening slightly the transition. Close to 
half-flux no commensurate state is available, thus phase 
correlations between network sites cannot be established, 
leading to a strong broadening of the transition. 

In fact, the T3 tiling geometry can be viewed as an 
ensemble of three coupled triangular sublattices, two 
formed by the 3-fold sites and another by the 6-fold sites. 
The singular properties of the T 3 spectrum e(/), at frus- 
tration / are simply revealed by the transformation: 

£ 2 (/)-6 = 2cos(7r/) e T (f T ) (2) 

that relates e(/) to the triangular lattice eigenvalues 
£ t(/t) at frustration fa — 3//2 At half-flux, due to 
cancellation of the cos(7r/) prefactor, all the energy lev- 
els collapse into two highly degenerate discrete levels at 
e = ±V6, forming flat, non-dispersive bands, in addition 
to the e = flat band. Due to the mapping of the TB 
problem onto the linearized GL approach, the superfluid 
velocity can be expressed in terms of the group veloc- 
ity of the band spectrum close to the ground state. In 
the context of a superconducting wave function, a non- 
dispersive state cannot carry phase information through 
the network, contrary to a Bloch state. Therefore, crit- 
ical current measurements give information on the net- 
work ability to sustain a supercurrent, i.e., both a finite 
order parameter and a finite superfluid velocity. 

The critical current was studied as a function of field 
from the dynamic resistance characteristics vs increasing 
dc bias current at temperatures close to T c (0). The used 
criteria was the threshold current for which the dynamic 
resistance exceeds 0.2% R n . Within the sensitivity lim- 
its of our measurements, it corresponds to the maximum 
current that the circuit is able to carry without dissi- 
pation. To avoid heating effects due to feeding a large 
current, we used sample Star 20 with 23 cells (20 fan) 
width. The critical current density per wire J C (T, f) is 
obtained from the network critical current divided by the 
number of parallel wires (25) and the wire cross section. 

Close to T c , we expect the critical current to follow a 
3/2 power law that generalizes the depairing current of 
a one-dimensional superconducting wire fl4|| to a super- 
conducting network p5|: 



J c (T,f) = J n C(f) 
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where J n is the zero field depairing current density at 
T = K. The field dependent coefficient C(f) is derived 
from the band curvature, d 2 e g /dk 2 close to the ground 
state, e g (f) by, 
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In Fig. H is displayed the field dependence of J c , at 
T = 0.96 T c (0) (1.185 K). Sharp peaks are obtained for 
the same frustrations as the downward cusps observed in 
the transition line. The remarkable finding is the total 
absence of peak in the critical current at the lowest order 
rational / = 1/2, exhibiting a clear minimum at this 
field. For all studied temperatures the critical current 
was always found to exhibit the lowest values at / = 1/2. 
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FIG. 3. Critical current density of Star 20 as a function of 
/ at T=1.185 K. The depression at half-flux is the signature 
of the non-dispersive state. Inset: Theoretical values at the 
same temperature. 



In the inset of Fig. [ 
obtained using (||) and 
T = 0.96 T c (0) and J n = 



is displayed the theoretical J c 
1) for / = p/q, q < 50. We used 
7x 10 6 A cm -2 , estimated from 



a 3/2 power law fit of the zero field J C (T, 0) data. The 
critical current consists of successive 5-functions at ratio- 
nal frustration / = p/q except at 1/2. The highest values 
of J c were obtained for / = 0, g, i and the symmetric val- 
ues. For these rational frustrations the spectrum is band- 
like and the group velocity is finite (Bloch states). The 
same applies for other regular periodic lattices such as the 
square lattice where a checkerboard commensurate state 
leads to an important peak at f=l/2 ]l5|] . However, in 
the T3 case at half-flux the group velocity is strictly zero 
due to the absence of dispersive states and the critical 
current vanishes. This situation is original for an infinite 
tiling and is due to the special T3 geometry. A similar 
effect is found when / approaches irrational frustration 
(for example, at small frustration / = l/q, with large 
q) : the bandwidth then becomes exponentially small and 
therefore the group velocity and the critical current are 
suppressed. 

The experimental data follow the same qualitative 
behavior vs field as the theoretical predictions. The 
3/2 power dependence of J c vs temperature at constant 
magnetic field was observed in the temperature range 
T c (f) - T < 20 mK. Namely, the field dependent coeffi- 
cient C(f) at half-flux is reduced to 17% of its zero field 
value. The C(f) depression, thus of J c , reflects the effect 
of the band structure on the superfluid velocity, and pro- 
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vides a strong evidence of the non-dispersive character 
of the state at / = 1/2, although the measured critical 
current does not vanish. 

One possible explanation for the incomplete suppres- 
sion of J c is the network finite size. A current carrying 
state (an edge state similar to surface superconductivity 
in finite type-II superconductors) exists along each edge 
of the finite network and is expected to lead to a non- 
zero supercurrent. A second possible origin for finite J c 
is the influence of the GL non-linear term which was ne- 
glected in Eq. (^|). Presumably, the non-linear terms in 
the GL formulation are responsible for degrading the fine 
features of the band structure and therefore give a finite 
critical current. The critical current observed in Fig. ^ 
at small frustrations, for example close to zero, may have 
the same origin. To go further, an exact solution of the 
non-linear GL equations would be needed. Nevertheless, 
as demonstrated by Abrikosov |jT6fl , a good physical in- 
sight of the superconducting properties can be obtained 
from the eigenstates of the linearized GL equation. 

This phenomenon suggests interesting properties of the 
vortex sublattice. In this context, the coupling between 
network sites can be expressed as a landscape of energy 
barriers against vortex motion. For example, at / = 1/3, 
a periodic vortex configuration can be easily constructed, 
matching perfectly the underlying lattice. This configu- 
ration is strongly pinned and very stable against driving 
currents, leading to a large critical current. The decou- 
pling of some network sites at / = 1/2 suggests that, in 
the absence of pinning, the vortex configuration will be 
highly disordered. Therefore, a significant dissipation is 
expected for small driving currents, as revealed on our 
experiments by the suppression of critical current and 
the anomalous transition broadening. These considera- 
tions are supported by preliminary experiments on vor- 
tex decoration which indicate a highly disordered vortex 
distribution at / = 1/2 and will be addressed elsewhere. 

More subtle is the commensurate state at / = 1/6, 
which corresponds to the fx — 1/4 state of the trian- 
gular lattice formed by the 6-fold sites (see Eq. ||). As 
shown in Ref. Ilffl, the uniformly frustrated XY model 
on a triangular lattice at fx — 1 /4 presents an accidental 
degeneracy of the ground state with zero energy domain 
walls which can weaken the global phase coherence. In 
our experiments we do observe a critical current peak at 
/ = 1/6, almost as large as at / = 1/3. The singular be- 
havior observed experimentally at / = 1/2 (fx = 3/4) is 
completely absent at / = 1/6 (fx = 1/4) and therefore, 
cannot be simply related to the triangular lattice prob- 
lem. Besides, it is not clear if the accidental degeneracy 
persists for a tight binding coupling. 

In summary, the anomalous transport behavior of the 
T3 superconducting networks at half-flux is consistent 
with the localization effect predicted in Ref. Q). The 
transition line is in excellent agreement with the related 
T3 ground state. The broad transition width at half- 



flux, indicates a strong enhancement of phase fluctua- 
tions which we assign to destructive quantum interfer- 
ences at this field. The reduction of the critical current 
at / = 1/2, which was never observed so far in periodic 
superconducting networks, illustrates the inability of the 
network to sustain a transport current. This behavior 
is the analog, in the superconductor case, of the metal- 
insulator transition predicted in Ref. 
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